CO VARIANT REPRESENTATIONS FOR MATRIX- VALUED TRANSFER 

OPERATORS 



DORIN ERVIN DUTKAY AND KJETIL R0YSLAND 

Abstract. Motivated by the multivariate wavelet theory, and by the spectral theory of transfer operators, 
we construct an abstract affine structure and a multiresolution associated to a matrix- valued weight. We 
describe the one-to-one correspondence between the commutant of this structure and the fixed points of the 
transfer operator. We show how the covariant representation can be realized on K n if the weight satisfies 
some low-pass condition. 
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1. Introduction 

Several themes are involved in this paper: the "covariant representations" in the title has been a central 
construct from the theory of operator algebras since the 1950s, and they have played a key role in numerous 
applications since. One of these more recent applications is to a class of wavelets called "frequency localized" 
wavelets. The "transfer operators" in the title refers to a construction with origins in probabilistic path 
models from physics and Ergodic theory. 

One of our aims here is to point out some connections between the two areas, and to show how operator 
algebraic ideas and representations throw new light on a classical theme. 

Since several ideas are involved, readers from one area may look for pointers to the other. 

We begin with a brief guide to the literature: One use of operator algebras (specifically, C*-algebras) is to 
the construction of representations. Initially [Mac78 , the focus was on groups, but the notion of covariance 
from physics (see e.g., |BR87I IRue04j ) suggested crossed products of groups with act by automorphisms 
on C*-algebras ([CMW84, Wil82 ). Since the pioneering paper by Stinespring (Sti55j , a preferred approach 
(e.g., |Arv 69 ) to constructing representations begins with a positive operator valued mapping, and it was 
Stinespring's insight that identified the correct "positivity" as complete positivity. But, at the same time, 
notions of positivity are central in a variety of probabilistic path models, beginning with Doeblin [Doe40] . see 
also |Coh81j . It is now also a key tool in ergodic theory, |Pet83| IWal82'J. As a result, Doeblin's operator has 
subsequently taken on a variety of other incarnations and it is currently known as "the transfer" operator, 
the Ruelle operator, or the Perron- Frobenius-Ruelle operator jB alOO] . The name Ruelle is from its use in 
statistical mechanics as pioneered by David Ruelle, see Rue04, BalOO . 
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Of a more recent vintage are applications to wavelets [Dau92] , i.e., special and computational bases in 
Hilbert space constructed from a class of unitary representation of certain discrete groups of affine transfor- 
mations. It was realized (e.g., |BJ02| ) that there are transfer operators Rw for wavelets, that the solution 
to a spectral problem for Rw yields wavelet representations; and moreover that these representations come 
with a useful covariance. Hence the circle closes with the positivity question from |Sti55j . albeit in a different 
guise. It is intriguin g that a variant of the operators Rw have now also found use in quantum error correction 



codes, see [CKZ061 lKrT06] . 



All the versions of transfer operators involve hierarchical processes with branching, and probabilities 
assigned by a weight function. Our paper focuses on Rw in the wavelet context; and we demonstrate that 
for many wavelets, the weight function must take the form of an operator transformation X i— » mXm* 
where m is a fixed matrix valued function and where m* denotes the adjoint operator, in this case transpose- 
conjugate. This matrix version of Rw is necessary for understanding wavelet constructions associated to 
wavelet sets |DL98j . and more generally to non-MRA wavelets, [BCM02] , (MRA stands for multiresolution 
analysis |Dau92j .) 

An orthogonal wavelet is a function ip € L 2 (R) such that {2^ 2 ip(2 : > ■ —k) \j,k 6 Z} is an orthonormal 
basis in L 2 (R). The main technique used to construct wavelets is by multiresolutions. For this one needs a 
low-pass filter mo, i.e., a 27r-periodic Lipschitz continuous function that satisfies the low-pass condition and 
the QMF condition: 

mo(O) = V2, 



mo 



(1) 



toq 
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Then, from the low-pass filter, one constructs the Fourier transform of the scaling function 



-r \ TT m ° (#) i 

71=1 ^ 



Finally, the wavelet is constructed from the scaling function using the formula 

= Tl 6 ' 1 ™° (1 + n ) * (1 ) ' (a;eM )- 

It turns out that the low-pass condition and the QMF condition, while necessary, are not always sufficient 
to obtain an orthonormal wavelet. One of the extra conditions on mo that guarantees the orthogonality of 
the wavelet was given by W. Lawton ; Law9T] . Here the transfer operator was introduced in the study of 
wavelets. 

In this context, the transfer operator is defined on 27r-periodic functions by: 

R m J(x) = i(K(|)| 2 /(f ) + K(^)l 2 /(^^)), (* e R). 

Lawton's condition states that the wavelet is orthogonal if and only if the only continuous functions h 
with R mo h — h are the constants. 

When this condition is not satisfied, the resulting wavelet still has an interesting property, namely it 
generates a Parseval frame. Having this, the theory of dilations of Parseval frames due to D. Han and D. 
Larson [HLOOj can be used. The wavelet Parseval frame is the projection of an orthonormal basis in a bigger 
space. The problem then was if this orthonormal basis has a similar wavelet structure, i.e., if it is generated 
by the application of unitary dilation and translation operators U and T that satisfy the commutation 
relation UTU~ X = T 2 (resembling the relation between the dilation and the translation operators in L 2 (R)). 
The answer to this question is positive and it is given in [BDP05j . It turns out that each fixed point of the 
transfer operator R mo will give rise to such a wavelet structure, and we call this a covariant representation. 
Putting together these representations, one obtains the covariant representation which has an orthonormal 
wavelet in a bigger Hilbert space whose projection onto L 2 {R) is the Parseval wavelet frame constructed 
in the classical wavelet theory. This is one of the wonderful uses of covariant representations. Since these 
orthonormal wavelet bases live in a bigger Hilbert space, i.e., one that contains L 2 (R) as a subspace, Han 
and Larson coined the term "super-wavelets" . Another interesting application of covariant representation 
is the computation of the peripheral spectrum for the transfer operator: the fixed points of R m „ are in 
one-to-one correspondence with the commutant of the covariant representation. Often this commutant can 
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be explicitly computed, hence the eigenspace of R mo can be obtained from that. The spectral properties of 
transfer operators play an important role in the ergodic analysis of discrete dynamical system (see [BalOO ) . 

Why matrix- valued transfer operators? It is known (see |Dau92p that not all wavelets in L 2 (R) can be 
constructed from a multiresolution. However there are generalizations of this that will do the job: each 
orthogonal wavelet can be constructed from a generalized multiresolution analysis, a notion introduced by 
Baggett et al. [BMM99J . The construction requires some matricial low-pass filters. Also multiwavelet theory, 
in which one uses more than one function to generate the basis, requires matricial filters. But many times 
the resulting wavelet is only a Parseval frame, not an orthogonal basis (see [BJMP05] ). 

There is no analogue for the Lawton condition in the case of matricial low-pass filters. The main reason for 
this is the impossibility to give a good generalization of the notion of zeros for the fixed points of the transfer 
operator, a notion which plays an essential role in the scalar case. We believe that the covariant representa- 
tions can provide a way around this, and we can analyze the orthogonality of wavelets and scaling functions 
constructed from matricial low-pass filters through a study of the associated covariant representations. 

More generally, a transfer operator, also called Ruelle operator, is associated to a finite-to-one continuous 
endomorphism on a compact metric space r : X — > X and a weight function W : X — > [0, oo), and it is 
defined by 

R w f{x)= ]T W(y)f(y) 

r(y)=x 

for functions / on X. 

Transfer operators have been extensively used in the analysis of discrete dynamical systems [BalOO] and 
in wavelet theory |BJ02| . 

In multivariate wavelet theory (see for example JS99J for details) one has an expansive n x n integer 
matrix A, i.e., all eigenvalues A have |A| > 1, and a multiresolution structure on i 2 (R"), i.e., a sequence of 
subspaces {Vj}j^z of L 2 (M. n ) such that 

(i) Vj C V j+1 , for aU j; 

(ii) UjVj is dense in L 2 (R n ); 

(iii) = {0}; 

(iv) / G Vj if and only if f((A T )~ 1 -) G Vj-x; 

(v) There exist ipi, (fd G Vo such that {<fk(- — j) \ k G {1, d},j G Z™} forms an orthonormal basis 
for Vq. 

The functions (pi,...,(pd are called scaling functions, and their Fourier transforms satisfy the following 
scaling equation: 

d 

<Pi&) = Q^mjiiA^x^iA-^x), (x G K", * G {1, d}), 
j'=i 

where q := | det A\ and rriji are some Z ra -periodic functions on M. n . 

The orthogonality of the translates of ipi implies the following QMF equation: 

- m*(y)m(y) = l, (x G M n /Z n ), 

^ Ay— x mod Z d 

where m is the d x d matrix (rriij)f - =1 . When the translates of the scaling functions are not necessarily 
orthogonal, one still obtains the following relation: if we denote by 

hijix) := &( x + k )^3^ x + ( x € K"/Z n ), 
fceZ" 

then the matrix h — (hij)f - =1 satisfies the following property: 

(1.1) Rh(x):=- V m*(y)h(y)m(y) = h(x), (x e M"/Z"), 

q * — ' 

Ay—x mod X n 

i.e., h is a fixed point for the matrix-valued transfer operator R. The fixed points of a transfer operator are 
also called harmonic functions for this operator. Thus the orthogonality properties of the scaling functions 
are directly related to the spectral properties of the transfer operator R. 

This motivates our study of the harmonic functions for a matrix-valued transfer operator. The one- 
dimensional case (numbers instead of matrices) was studied in |BJ02[ IDut04b[ IDut04a| . These results were 
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then extended in [DJ06a, DJ06b|, by replacing the map x i— > Axvaodl/ 1 on the torus T"/Z n , by some 
expansive endomorphism r on a metric space. 

Here we are interested in the case when the weights defining the transfer operator are matrices, just as in 
equation (II. 1|) . We keep a higher level of generality because of possible applications outside wavelet theory, 
in areas such as dynamical systems or fractals (see |DJ06al rDJ06bj ). However, for clarity, the reader should 
always have the main example in mind, where r : x i— > Ax mod Z™ on the torus T™ . 

In |DJ0 6b] it was shown that, whenever a pair (m, h) is given, with h > and Rh = h, one can construct 
a covariant structure on some Hilbert space H, i.e., a unitary U , a representation 7r of continuous functions 
on X, and some scaling functions ip±, . . . , ipd € H such that 

U*(f)ir=n(for), (feC(X)), 

</>,•) - / (/eCpO), 

(1.2) t/^i = ^7r(m^)^-, (*G{l,...,d}), 

(1.3) {{/- n 7r(/)^ |n€Z,ie{l,...,d},/€ C(X)} is dense in if, 

where /i is a strongly invariant measure on X (see (12.11) ). 

Moreover, this construction is unique up to isomorphism. Thus, from the "filter" m and the "harmonic 
function" h, one can construct a multiresolution structure, similar to the one used in wavelet theory (see 
BJ02 ). There are several uses for such a structure: one can analyze the peripheral spectrum of the transfer 
operator Dut04bJ, construct super- wavelet bases on spaces bigger than L 2 (R) or on some fractal spaces 
BDP05, DJ06c], or use the rich algebraic and analytic structure of this multiresolution to perform compu- 
tations needed in the harmonic analysis of fractal measures |DJ06a| . 

This covariant representation was analyzed in more detail in the one-dimensional case d = 1 in [DJ06b, 
D J06a] . The main tool used there was the introduction of some random- walk measures P x following an idea 
of Conze and Raugi |CR90j . 

In our case m and h are matrix- valued. We will use the language of vector bundles because we are interested 
also in continuous harmonic functions and projective multiresolution analyses (see [DR061 lPR .04 _). We will 
describe the covariant representation in the multivariable case using some positive-matrix valued measures 
(see Theorem [321 Propositions 15.11 and I5.3[) . thus extending previous results from [DJ06bJ. 

In Section [2] we introduce the setup and the main definitions and assumptions on the transfer operator 
R associated to a matrix- valued weight h. Then in Section [3] we show how the covariant representation 
can be constructed on the solenoid of r. In Section 0] we show that the operators that commute with this 
covariant representation are in one-to-one correspondence with the harmonic functions of R (Theorem 14. ip . 
Then the set of harmonic functions inherits a C*-algebra structure from the commutant, and a more intrinsic 
description of the multiplication is given in Theorem 14.31 and Remark 14.41 

The covariant structure can be described in terms of some operator- valued measures P x . We do this in 
Section [5] With the aid of these measures we can give another form to the correspondence between harmonic 
maps and operators in the commutant, which can be identified now with cocycles (Proposition 15 .61 Theorem 

In the case when the filter m satisfies a low-pass condition (in this case, the £J(Z)-condition as defined in 
JS99]), we can realize the covariant representation in the more familiar environment of i 2 (K"). The space 
W 1 has a natural embedding in the solenoid, the measures P x are atomic and, the measure of the atoms are 
directly related to the solutions of the refinement equation, i.e., the scaling functions. Our results extend 
some ideas from |Gun00i IDGHOOi IJor06j to the matricial case. 

2. Definitions and preliminaries 

The dynamical system. Let X be a compact Hausdorff space with a surjective and finite to one 
continuous map r : X — > X . Moreover, let /i be a regular measure on /j, that is strongly r- invariant, i.e. 
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for every f £ C(X). 

Main example. The main example we have in mind, is the following: Let G C l n be a discrete subgroup 
such that R"/G is compact, i.e., G is a full-rank lattice and W l /G ~ T n =: X. Moreover, let A £ GL(R n ) 
be strictly expansive and such that AG C G. Let p : E™ — * X denote the quotient map and define a map 
r : X — > X as r(p(x)) = p(Aa;). This is a | det ( A) | -folded normal covering map. Finally, we let /i be the 
Haar measure on T™. 

The transfer operator. Let p : £ — > X be a <i-dimensional complex vector bundle over X with a 
Hermitian metric on £, i.e., a continuous map (•, ■) : £ x £ — > C that restricts to an inner-product on each 
fiber. Such a map always exist when X is compact, [Ati89 , 1.3.1]. Let S be the set of continuous sections in 
£. Then 

•si,s 2 ^ / (si(x),s 2 (x))dfj,(x), 

JX 

defines an inner-product on S. Let K, denote the Hilbert space obtained by the completion with respect to 
the corresponding norm. Then JC is the Hilbert space of L 2 -sections in £ with respect to the measure /z. 

The C(A)-product on the sections in £ gives a representation of C{X) on K, by pointwise multiplication, 
k : C(X) — > B(JC) such that (n(f)s)(x) = f(x)s(x) for every x £ X. Every bundlemap on £ commutes with 
this representation. 

Remark 2.1. Let us consider the von Neumann algebra k(C(X))" generated by this representation. Since 
the bundle is locally trivial around any point in X, we can take U\, . . . , U r open sets in X, with bundle 
isomorphisms <j>j : £\u, — > Uj x C d . 

If V £ k(C(X))" there exists V(x) £ C such that cc i— > (j)j(x)~ 1 V(x)<pj(x) is an L°°(Uj, function and 
0j(Vs(a;)) = 4>j(V(x)s(x)), /x-a.e. 

If y G k(C(X))' there exists V^(ir) € Endc(£|a;) such that the map x i— > :) (a;)~ 1 V r (x)(/) :) (x) is contained in 
M d (C) ® L^E^i/lo,) and ^-(^s)(a;)) = <f> j {V(x)a(x)) fx-a,.e. 

Let r*£ denote the pull-back of £ along r (see [Ati89, 1.1]). The space of sections in r*£ is endowed 
with the pull-back Hermitian metric from £. Let K, denote the L 2 -sections in this bundle with respect to 
this metric and /i. Moreover, let k : C(X) ~ * B{JC) denote the corresponding representation of C(X) by 
pointwise multiplication. 

The weight (or the fitter if we use wavelet terminology) that is used to define the transfer operator is 
in our case an operator m 6 B(K,,K,) such that mk(f) = n(f)m for every / £ C(X). As in Remark 12.11 
m is a pointwise multiplication by a linear map between the fibers £ rx and i.e., there exist a unique 
m(x) € Home (£ |ra:, £U) sucn that (f>j(m(x,v)) — <fij(m(x)v) for every a; € J7j and u € ( rl . 

Definition 2.2. Let = k(C(X))' . As we have seen before, the space M. consists in bounded measurable 
bundle maps on £. We define the transfer operator on M by 

(Rf)(x) = 1 m *(v)f(y) m (y)> 

TT \ J r y — X 

Assumptions on the filter m. Throughout the paper we will assume that m satisfies the following 
conditions: 

(i) m 6 i?(/C,/C) is injective 

(ii) sup fc ||i? fc || < oo 

(iii) There exists an h G End(£) such that Rh — h and h > (here we mean h non- negative as an 
operator). 

Remark 2.3. The first condition ||]J implies that m(x) is invertible /i-a.e. 

Note that if there exists an h G M. such that Rh = h and h > cl for some c > (here 1 is the identity 
bundle map on £), then ((n]) holds automatically, sec DR06J. 

Definition 2.4. A function f £ M. such that Rf = / is called harmonic with respect to the transfer 
operator R. We denote by the set of all bounded harmonic functions. 
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3. The COVARIANT REPRESENTATION associated to a matricial weight and a harmonic bundle 

MAP 

We mentioned in the introduction that for some choices of the filter mo such as the "stretched Haar filter" 
(see |Dau92j ) the scaling function and wavelet constructed in L 2 (M) are not orthogonal, i.e., they do not have 
orthogonal translates. Thus classical wavelet theory breaks down if we want to have orthogonal solutions for 
the prescribed scaling equation. Even deeper problems appear when the filter does not satisfy the low-pass 
condition, since the infinite product defining the scaling function from too is 0. Fractal spaces may occur 
(sec [D J06c ) . We will show that we always have a solution to a given scaling equation in some Hilbert space 
with a covariant representation. The Hilbert space depends sensitively on the filter m and the correlation 
function h. 

We construct now the affine structure associated to the matricial filter to and the positive harmonic 
function h. The covariant representation is a Hilbert space endowed with an affine structure given by a 
unitary U, which takes the place of the dilation operator, and a representation n of C(X), which takes the 
place of the representation generated by translations. In this Hilbert space one has several vectors (pi, . . . , ip n 
that satisfy the scaling equation with the prescribed filter to as in (|1.2p . and having correlation function h 
as in p.3[) . 

The ground space. Consider the projective system 

r v 

X* X* X* 



and let X^, be its projective limit, with projections 9 n : — > X, (n > 0), and a homeomorphism 
f : X x — > Xoo such that the following diagrams are commutative for all n > 0: 




We have the following identification of X^ 

Xoo = {{x$, xi, X2, ■ ■ ■ ) € X X X X . . . \rxj+\ = Xj, for all j > 0} 

9k(xo,xi,..-) =x k , f{xQ,x\,...) = (rxo,XQ,xi,...). 
The space X^ is a compact Hausdorff space with the topology generated by the inverse images of the open 
sets in X with respect to the maps 6 n , n > 0. 

The Hilbert space. Let Sk denote the vector space of continuous sections in 6q£ that depend only on 
the first k + 1 coordinates Xo, x\, . . . , Xk- We see that Sk is formed by s o 0^ for sections in £, s : X — > £, i.e., 
maps of the form 

x h-> (r~ k (x), s o 9 k (x)), (x e X). 
We define the inner-product (•, : Sk x Sk — * C by 

(fo9 k ,go9 k )k ■= j x #r l k{x) E {m ik Hy)f(y),h(y)m^(y)g(y))d^x), 

r k y—x 

where 

TO,' fe '(a;) := m(x)m{rx) . . . m(r k ~ 1 x). 
We have the following compatibility relation between these inner-products: 

{f°9i,godi)k+i =(f°r k o9 k +i, g ° r k o k +i) k+l 



#r -(l +l){x) E {™ {k+l \y)f{r k y)MyW k+l) {y)g{r k y ))d^x) 



v ' r K ~ l ~ L y—x 



l x #^T) E {m {l \y)f{y),R k (h)(yW l) {y)g{y))dn{x) 



r L y—x 



=(/ o9 h go9()i (since R k h = h). 
This shows that the restriction (•, ^k+ilsixSi — {■, •);• 
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We obtain a sesquilinear form (• , •) on U k S k , i.e., a (possibly degenerate) inner-product. Let TL denote 
the Hubert space completion with respect to this inner-product. 

An application of the Stone- Weierstrass theorem gives that U k S k is dense in the space of continuous 
sections in 9q£. Moreover, the continuous sections in are contained in TL by the following argument. 

Lemma 3.1. If {sk}keN is a sequence of sections in £ such that {s/c°#fc}fc>o converges uniformly to a section 
s : Xoo — > #q£ then (s, s) = lim fc (s fc ° S k ,s k oO k ) k . 

Proof. 

\\s k+ i o e k+ i - s k o 9 k \\ 2 

= / #r -(*U( x ) £ {Sk+l{y) - Sk ° r l (y)M k+l) *{y)h{yW k+l \y){ Sk+l (y) - s k o r l {y)))d^x) 

r k + ly — x 

<\\ Sk+l o e k+l - Sk o OuWl J #r _ (fc 1 +0(a;) E \\m {k+l) *{y)HyW h+l ){ y )\W{x) 

r k+ly — j. 

<\\s k+ i o 9 k+l - s k o ^ll^rfll^+^l^, 

and the last inequality follows since £\ \\HIH.W < £\ Tr(P*P,) = TrQ^P*^) < d\\ J2i H t H i\\ for ever y 
finite sequence of matrices Hi. □ 

Next we define the covariant representation and the multiresolution on the Hilbert space TL. 

Theorem 3.2. Let P k denote the orthogonal projection in B(TL) onto the closed subspace TL k generated by 
S k . Define 

7T : C(X) — > B(TL), (n(f)s)(x) = f o o (x)a{x), (/ G C(X), s G TL,x G X^); 

U G B(TL), (Us)(x) =mo6 {x)sor(x), (s€TL,x€X). 
The following relations hold: 

(i) U is a unitary such that Un(f)U* — 7r(/ o ?-) for all f G C(X). 

(ii) P G n(C(X))'. 

(iii) P k < Pfc+i for every k G N. 

(iv) limj. P k s = s for every s G TL. 

(v) UP k+1 U* = P k for every k G N. 



Proof. dTTJ) and (JTTTJ) follow directly from the definition. (fry)) follows from the fact that U k S k is dense in the 
sections of 9q£ . 

We must prove (0). First, let f,g G 5, then 

(m o 6> / o9 k of,mo 9 g o 6> fc of) = ((m o r fc_1 /) ° <9fc x , (m o r k ~ 1 g) o (9 fc _i) 

„ J Iw ; E (m^- 1 \y)m(r k - 1 y)f(y),h{y)m^- 1 '>(y)m(r k - 1 y)g{y))d t x{x) 



1 ^ 'm^ (y)f(y), h{y)m^ (y)g(y))dfi(x) 



#r-(*-V(x) f 

r fc a.— . 

=(foe k ,go9 k ) 

This shows that [/ is isometric on UfcSfc. 

To see that J7 is surjective on TL, we will show that UTL k +i = 7Yfc, (with 7Yfc := Pfc7Y) which by ((TvJ) 
implies that [7 is surjective. If s G TL k +i, then J7s — mo 9qs org 7ifc. If we can show that {(m o r k s) o 
k \s measurable section in £} is dense in TL k , we are done. Recall that m(x) is invertible /i-a.e. 

Take s E S. Suppose, is a decreasing family of open sets such that {x G X| det to or fe (x) = 0} C Vi 

for every / G N and lim;/i(V;) = 0. (Note also that by (|2.ip the measure /i is invariant under r, i.e., 
fi(r (E)) = jii(-E) for all measurable sets E). For every Z G N, there exists a measurable section si in £, 
such that s; is on VJ, and m o r k (x)si(x) = s(x) for .t G X \ Vi. Then lim; to o r fc s; = s, /i-a.e. and by the 
dominated convergence theorem, lim; Usi o 0^+1 = s o 9 k , in 7i, i.e., s o 9 k is in the closure of U S k . 

Finally (jvj) follows from the relation UTL k +\ — TL k and ([!]). □ 
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Definition 3.3. We denote by C*(X, r, m, h) the C*-algebra generated by U and 7r(/), / G C(X) 1 and we 
call it the covariant representation associated to m and h. 

4. Harmonic and measurable bundle maps 

In this section we establish the one-to-one correspondence between harmonic functions and operators in 
the commutant C*{X,r,m,h)' . We begin with a result on self-adjoint elements, and we require that the 
harmonic function ho be dominated by h, in the sense that \ho\ < ch for some constant c > 0. In the case 
when h is bounded away from 0, i.e., h > CI for some C > 0, then this domination condition is automatically 
satisfied and we can extend the correspondence to non-self-adjoint elements (Corollarv l4.2p . 

Theorem 4.1. (i) [Harmonic maps to operators in the commutant] If ho G A4, h^ = ho, Rho = ho 

and there exists a positive number c > such that \ho\ < ch, there exists a unique self adjoint operator A in 
the commutant C*(X,r,m,h)' such that 

(f o O , Ag o 9 ) =J(f, h o9) , (f, 9 e S). 

We denote this operator by Ah Q ■ Moreover in this case 

(fo9 k ,Ago9 k ) = [ 1 T (m^(y)f(y),h (y)m^(y)g(y))n(dx), 
Jx #r k {x) 

r K y—x 

for every k € N and f,g € S. 

(ii) [Operators in the commutant to harmonic maps] Conversely, define the operator T G B(PoH, JC) 
such that Tso0 {) = hs for every s € S. If A is a self-adjoint operator in the commutant C*(X,r,m,h) then 
hA '■= TPqAPoT* G A4 is a harmonic function, i.e., PJia — hA, and \Ha\ < ch for some constant c > 0. 
Moreover, the correspondences described in (i) and (ii) are inverses to each other, i.e., 

h AhQ = ho, A hA = A. 
Proof, (i) Let R>k denote the sesquilinear form on Sk, given by 

B k (foe k ,go9 k )= [ - 1 V ( m M(y)f(y),h (y)mM(y)g(y))n(dx). 
Jx # r W k 

Using the computation for (•, -} k , we see that B k \s k _ 1 xs k -i — B k -i- 

By the boundedness assumption on ho, we have —ch < ho < ch. This shows that 

-c (/ o 9 k , / o 9 k ) < B k (f o9 k ,fo 9 k ) <c(fo9 k ,fo9 k ). 

Since ho is self-adjoint this implies that we obtain a bounded sesquilinear map B on Ti that restricts to B k 
on S k . Let A denote the bounded operator on H such that B(f, g) = (f, Ag) for every f,g&H. 

The computation that showed that U was isometric with respect to (•, •) applies here too, and it shows 
that B(f, g) = B(Uf, Ug) for every f,g€H. Now 

(/, UAg) = (U*f, Ag) - B(U*f, g) - B(f, Ug) = (/, AUg), 

i.e., AU = UA. Moreover, by a direct computation, we have B(f, ir(a)g) = B(n(a*)f, g) for every a € C(X), 
f, g G Ti.. Using this, we obtain 

(f,Air(a)g) = B(f,ir(a)g) = B(Tr(a*)f,9) = (ir{a*)f,Ag), 

i.e., ir(a)A = Air(a) for every a G C(X). Finally, since U is unitary, U* A = U~ 1 A = AU^ 1 = AU* and A 
commutes with every operator in C*(X, r, m, h). 

The uniqueness follows from the fact that U k S k is dense in Ti. 

(ii) First, a simple computation shows that T* s = s o 9o for s 6 K. Then 

/ (si(x), (TaT*)(x)s 2 (x))dfi(x) = (s 1 o 9 ,as 2 o 9 ). 
Jx 

for every a G n(C(X))' D B(P Ti). 

If A commutes with 7r(/) for all / G C(X) then, since Po commutes also with n(f), it follows that PoAPo 
commutes with tt(/) so TP AP T* commutes with «(/) for all / G C{X). Therefore TPqAPoT* G k(C(X))', 
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and we can use Remark |2. II Moreover, if A € C*(X, r, m, h)', it commutes also with U , and then, using the 
strong invariance of the measure \x, 

/ ( Sl ,R(TP AP T*)s 2 )d»= f -^J2( m (y) s dry),(TPoAP T*)(y)m(y) S2 (ry))d^x) 
Jx Jx W x ry=x 



= / (m(x)si(rx), (TP Q AP T*){x)m{x)s 2 {rx)}du(x) = (Usi o 9 Q ,P Q AP Q Us 2 o 9 ) 
Jx 

= (Us 1 o9o,AUs 2 o6o) = (s 1 o0 Ot As 2 o0 o ) = (si o o> P o AP o s 2 o O ) 

= [ ( Sl ,TP a AP T*s 2 )dfi, 
Jx 



i.e., TP AP T* is a fixed point for R. 

Since A is bounded and self-adjoint, it follows that 



(si(x) , h A {x)si(x)) dfi(x)\ = | (si o O , Asi o 9 ) 



x 



< \\A\\ (s a o 6»o , si o 6» ) = ||A|| / (si(a;) , /i(x)si(x)) d/x(x), 

Jx 

which implies that — \\A\\h < Ha < ||^.||^» f° r A*-a.e. a: € A". 

To prove that the correspondences are inverses to each other we compute: 

/ ( Sl ,TP AP T* S2 )d^= ( Sl o0 ,P AP s 2 o6 ) = ( Sl o9 , As 2 o6 ) = [ (s 1 ,h s 2 )dn, 
Jx Jx 

for every Si,s 2 € E, i.e., h = TP AP T* and A = A ho . 

□ 

Corollary 4.2. Assume in addition that h > cl for some constant c > 0. Let T G B(PqH.,]C) such 
that Ts o 9q — hs for every s G S . The map A i— > TPo^-Po? 1 * yields a bijection between the commutant 
C*(X, r, m, h) 1 and the bounded harmonic functions $j := {/ G M \ Rf = /}. 

Proof. Since h > cl, the boundedness condition \ho\ < const h is equivalent to the essential boundedness of 
ho, therefore it is automatically satisfied. Theorem 14. II gives us the bijection between self-adjoint elements. 
But every element A in the commutant can be written as a linear combination of self-adjoint elements: 
A = (A + A*)/2 + i(A — A*)/2i. Similarly for the harmonic functions. This gives the bijection. □ 

Corollary 14.21 shows that the set of bounded harmonic maps has also multiplicative structure, the one 
induced from multiplication of operators via the given bijection. The following theorem gives an alternative, 
more intrinsic description of this multiplicative structure. It is also a generalization of a result from |Rau94j . 

Theorem 4.3. Suppose h > cl [i-a.e., for some constant c > 0. If h±,h 2 G f) then 

hi * h 2 (x) := lim R k (h 1 h~ 1 h 2 )(x) 

k 

exists fi a. e and hi * h 2 G io. 

The map A i— * TPqAPqT* defines a ^-isomorphism from C*(X,r,m,h)' to i5, when S) is equipped with 
the product * and the ordinary involution. 

Proof. We define R(b) = h- 1 / 2 R(h 1/2 bh 1 / 2 )h- 1 / 2 . Now R is completely positive and unital (see |DR06j ). 
The Kadison-Schwarz inequality Bla06, II. 6. 8. 14] for completely positive and unital maps yields 

RibYh-'Rib) = h^Rih-^bh-^YRih-^bh- 1 / 2 )^/ 2 < h^Rih-^h^bh- 1 / 2 )^/ 2 = R{b*h- l b) 

This implies that R^h^hi) > R(h*)h~ 1 R(h i ) = h^h' 1 ^, i G {1,2}, so R k (h*h- 1 h l ) is a positive 
and increasing sequence in M. Since sup fc ||i? fe || < oo, the sequence is uniformly bounded, so R k (h*h~ 1 h i ) 
converges pointwise /z-a.e. to some element hi * hi in M.. And it is easy to see that Rihi * hi) = hi * hi. 

Note that hih~ 1 h 2 = jT,l=o i ~ k ( h i + ^M^^Oi +i k h 2 ). It follows that lim fe R k (h 1 h~ 1 h 2 )(x) exists 
/z-a.e. for arbitrary h±,h 2 G Sj. 

If Ai,A 2 G C*(X, r, m, h)' , hi := TP AiP T*, h 3 := TP a Ai A 2 P T*, then, since 

(si o 9 , A t s 2 o 9 ) = / (s 1 (x) , hi{x)s 2 (x)) du{x) = (s 1 o 9 , (h^ 1 ^) o 9 s 2 o 9 ) , 
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we see that PoAiPg is an operator of multiplication by hih 1 o O . 
Let a := hih~ l h 2 — h 3 . Then 

|| (P fe Ax P k A 2 P k - P k AiA 2 P k )f o O || 2 = \\{P A X P A 2 P - P A x A 2 P )U k f o O || 2 

= f (h- 1 (h 1 h- 1 h 2 -h 3 )m {k) for k ,hh- 1 (h 1 h- 1 h 2 -h 3 )m {k) for k )dn 
Jx 

= [ (f,R k (a*h- 1 a)f)d^ 
Jx 

since U is unitary, P k = U PqU , and /i is strongly invariant. 

Since P k converges to the identity strongly as k increases, we obtain that lim^ J x (f,R k (a*h~ 1 a)f)d/i = 
for every / G K, But this implies that there exists a subsequence such that linx,- R kj (a*h~ 1 a)(x) = /z-a.e. 

The Kadison-Schwarz inequality implies that R k (a)*h~ 1 R k (a) < R k (a*h~ 1 a). Then R k j (a) converges 
fi-a.e. to 0, so R kj (hih" 1 ^) converges /z-a.e. to h 3 . We already know that R (h x h h 2 ) converges fi-a..e. 
to hi * h 2 , so hi * h 2 (x) — h 3 (x), /i-a.e. 

□ 

Remark 4.4. We proved in UDR06J, that in the case when r is obtained by applying a certain covering 
projection to an expansive automorphism, and m and h are Lipschitz, with h > cl, then there is a C*- 
algebra structure on the continuous harmonic functions. The multiplication is constructed as follows: due 
to the quasi-compactness of the transfer operator R (restricted to Lipschitz functions), the uniform limit 

-. n— 1 

Ti(/):= km -Vff'/, 

n-^oo 77, ' ^ 

j=o 

exists for every continuous / and defines a continuous harmonic function. The product of two continuous 
harmonic functions is defined by (hi,h 2 ) i— > T^hih^ 1 ^). Theorem 14.31 shows then that this product 
coincides with hi * h 2 . In particular, if hi,h 2 are continuous harmonic functions, and A^.A^ are the 
associated operators in the commutant, then Ah 1 Ah 2 — Ah 1 *h 2 is a l so associated to a continuous harmonic 
function hi * h 2 = Ti^hih -1 ^). 

The next corollary shows how the covariant representation can be decomposed using projections in the 
algebra Sj. 

Corollary 4.5. Suppose h\,...,hi form a family of mutually orthogonal projections in fj, i.e., 

hi * hi = h* = hi, 1 < i < I, hi * hj = 8ij, 1 < i,j < I. 

Let pi, . . . ,pi denote the corresponding orthogonal projections in C*(X,r,m,h)' . Moreover, let (Wj, (•,•)«) 
denote the Hilbert space, associated to hi, in Theorem \ 3.'2\ with the representation TTi : C{X) — > B(7ii), 
unitary Ui and increasing family of projections {P k }k- 

Then p{Yi — Hi, (piSi,pi.s 2 ) — (si,s 2 )i for every Si,s 2 sections in 9q^. We obtain an isometry 

J : ®- =1 Wi -» H, J(si, ...,st)= pisi + . . .pisi, 

such that 

(i) n(a)J = J(® l i=1 n(a)) 

(ii) UJ = J(® l i=1 Ui) 

If^2ii—i hi = h then J is a unitary. 

5. COCYCLES 

In this section we give an alternative description of the Hilbert space Ti and of the operators in the 
commutant in terms of some matrix valued measures P x on the solenoid Xoo. 

In the scalar case, the measures P x are random walk measures with variable coefficients. Each point in 
the solenoid Xrx, can be constructed in the following way. Pick a point xo G X. Then, since rxi — xq, one 
has to make a choice xi out of the finitely many roots in r~ l x$. The QMF equation amounts to 

£ K( y )| 2 = i. 
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Therefore \mo(y)\ 2 / #r _1 (y) can be interpreted as the probability of transition from x to its root y. At the 
next step one makes a transition from x\ to its root X2 with probability given by \mo(x 2 )\ 2 / z fl z r~ 1 xi. And 
so on. The measure P x is the path measure obtained in this way. 

In our matricial case, the measures P x will be operator valued, but the idea stays the same. However 
some complications appear because we are dealing with a non-commutative situation. 

We will also see that the operators in the commutant of the covariant representation are in fact mul- 
tiplication operators by matrix valued functions. This will enable us to give a more concrete form of the 
correspondence between harmonic maps and the operators in the commutant. The result resembles the 
Poisson-Fatou-Privalov theorems in harmonic analysis: the harmonic map is the integral on the boundary 
of the operator in the commutant, and the operator is a radial limit of the harmonic map. In our case the 
boundary is the solenoid X^. 

First, by applying the Trace we will convert the matrix measures into a scalar measure and a matrix 
valued Radon-Nykodim derivative. 

Let Ck C C{X 00 ) be the set of continuous functions on Xoo that only depend on the fc + 1 first coordinates. 
Then 

f°0k" I Tv(R k (hf)(x))dn(x) 
Jx 

defines a positive linear functional <E such that the compatibility condition Tk+i\c k = Tk holds. 

Since sup fe ||-R fc | < oo and UfcCfc is dense in C(X ao ), there exists a positive functional r € C(X oc )* such 
that r\c k — Tfc for every fc. Let (i be the measure on Xoo, provided by the Riezs representation theorem such 
that t(/) = J x fdfi 

Proposition 5.1. (i) There exist a positive A(x) <G Endc{0^£\ x ), for fi-a.e. x 6 Xoo such that ||A(-)|| € 
L 00 (X 00 , jl) and 

(sx(x), A(x)s 2 (x))dp, = (si, s 2 ), 



for every pair of sections si, s 2 ■ Xoo — > #q£. 

(ii) If A is an operator in the commutant C*(X,r,m,h)' then for fi-a.e. x € Xoo, there exists A(x) £ 
Endc(9Q^\ x ), such that As (x) — A(x)s(x), fi-a.e., and 

(5.1) A(x)m o 6 (x)A(f(x)) = A(x)A(x)m o d Q (x), jl-a.e. 

Also, conversely, any such esentially bounded matrix-valued function x ^— > A(x), satisfying (I5.1[) . defines an 
operator in the commutant C*{X, r, m, h)' . 

Proof. Denote L 2 (9q£, fi) denote the L 2 -sections in the bundle 9q£ with respect to the measure fi. Note that 
if s is a section in £, then 

(so6 k ,so6 k ) = f V —±—{ m W(y)s(y),h(y)mW(y)s(y))dt, 

r K y—x 

E 4i l k, d m(ky(y)h(y)mW(y)\\(s(y),s(y))drix) 
x , #r L {x) 

E ^^^{rn {k) *{y)h{y) m ^\y)){s{y),s{y))d^x) 



(s o 6 kl s o 9 k )dfjL, 



x, 



so there exists a unique and positive bounded operator on L 2 (6*q£,/2) such that J x (si(x), (As 2 )(x))dfi = 
(si, s 2 ), for every pair of sections at, S2 ■ X30 — > #q£. 

The C(X 00 ) module structure on the section in 9q£ gives us representations of C(X oc ) on Ti and L 2 (9q£, fi) 
as follows: Every / S C(X oc ) gives a multiplication operator Mf € B(Tt) such that (Mfs)(x) = f(x)s(x). 
We see that Afjf = Mj as operators on both H and L 2 {9^, /}), so 

{at{x),(AM f 8 2 ){x))d[x= (at,M f a 2 ) = (Mj Sl ,s 2 ) 
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(M J s 1 (x),As 2 (x))dfl = / (s 1 (x),M f As 2 (x))dfl, 

A^oo J Aqo 

and therefore MfA — AMf. This implies that As(x) = A(x)s(x) for a A(x) S End(#Q£U) as described, for 
/t-a.e. a; € X. 

(ii) First note that C* (X,r, R, h)' C C(Xoo)'. This is because U^ 1 Tr(f)U = Mf g l: so if an operator 
commutes with U and 7r then it must commute with all multiplications by functions which depend only on 
finitely many coordinates. But these functions are dense in C{X oa ) and we obtain the inclusion. 

From this, we see that there exist A(x) £ #o£U sucn that As(x) — A(x)s(x), /2-a.e. Moreover, UA = AU 
implies the relation (|5.1[) . The converse follows by a computation to prove that A defined from x i— ► A{x) 
commutes with U and 7r(/) for all / £ C(X). □ 

Definition 5.2. Because of the relation (|5.ip . if A is an operator in the commutant C*(X, r, m, h)' then we 
call A a cocycle. 

In the next proposition we will define the positive-matrix-valued measures P x which can be used to 
represent the inner-product on the Hilbert space Tl. 

Proposition 5.3. For each x £ X, let Q x — {y £ X 00 \8o(y) — x}. There exists a positive operator-valued 
measure P x from the Borel sets in Vl x to Endc(£,\ x ) such that J„ / o 9kdP x = R k {fh){x), for every bounded 
measurable function f on X , and 



[si(y),dP x (y)s 2 {y))dfj,(x) = (si,s 2 ), 

a Jn^ 

for every pair of sections si, s 2 ■ X x — * $o£- 

Proof. Note that Si(y) £ 8q£\ x for every y £ £l x . Let Ck, x C C^^) denote the functions that only depend 
on the first k + 1 variables. Moreover, define : C(Vl x ) — > End^lx) by 

o- x k (fo6 k )=R k (fh)(x) 

We see that crjj; defines a positive bounded operator from C kyX to M<j(C). Moreover the compatibility condition 
a k+i\c k — °~ k holds. And, since sup fe ||i? fe || < 00 and \J k C k)X C C(fl x ) is dense, we see that there exists a 
positive and bounded linear map P x : C(Q X ) — > End(^| 2; ) for almost every x G X, such that cr^lc x = o%- 
A matrix computation implies that 

/ / (si o 9 k (y),dP x (y)s 2 o 9 k (y))dfi(x) = (si o 9 k , s 2 ° 6 k ), 

for every pair of sections si,s 2 : X — > £. The last claim follows from this by a density argument. □ 

The next Lemma can be obtained by direct computation. 
Lemma 5.4. For x £ X, let fi x be the measure on Q x defined by: 

f o 9 k dfi x = Ti(R k (fh)(x)) = Tr( f fo6 k dP x ). 
Then, for all bounded measurable functions on X^ 



fdfi= / f dfi x dfJ,(x). 
Aoo J x in, 

If Cwi,...,a)„ is the set of points (x, x\, x 2 , . . . ) £ Cl x such that x\ = uji, . . . , x n — uj n , then 

P*{C Uu ..^ n ) = #r _ 1 . MC^...,uJ = #r l n{x) Tr ((mWW"))K)) . 

Having defined the matrix valued measures P x , the correspondence between cocyles and harmonic func- 
tions in Theorem 14. If ii) can be given now in terms of a matrix valued conditional expectation: 

Proposition 5.5. Let A £ C*(X,r,m,h)' and A(x) £ End(9^\ x ) such that As(x) — A(x)s{x), fi-a.e., for 
every section s. We have the following identity: 

TP a AP Q T*(x) = I dP x (y)A(y), 
i.e., x 1— > J n dP x (y)A(y) is a fixed point for R. 
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Proof. If si, S2 ■ X — ► £ are sections, then 

(si{x),{TP Q AP T*)(x)s 2 (x))dfi(x) = (sio8o,P APoS2o8 ) 



x 



(sio9 ,As 2 o9 ) = I (si(x), I dP x {y)A(y)s 2 (x))dfi(x). 

□ 



For the inverse correspondence in Theorem I4.1f i), from harmonic functions to cocycles, we have the 
following result: 

Theorem 5.6. Assume in addition that h > cl, fi-a.e., for some constant c > 0. Let ho be a bounded 
harmonic function and let A be the corresponding cocycle as in Theorem \4- l^ i) and Provosition \5. 11 Then 

(mW*h mW) o0 k 
urn —j — = AA, 

fe^oo Tr((m( fc ' hm^ k i) o 9k) 

pointwise fi-a.e. 

Proof. For / bounded measurable function on let Ek(f) denote the conditional expectation onto the 
functions that depend only on the first k+1 coordinates, with respect to the measure fx. If F is a matrix-valued 
function on X^, then E k (F) is the matrix- valued function obtained by applying Ek to each component. 

Let Ak := (m^ h^homP^) o 9 k - (Recall that m(x) is invertible for /z-a.e. x G X). Then 

(/ o e k , A kg o e k ) = (/ o e k , a 9 e k ) = {fo9 k , p k AP k9 e k ) . 

Take f k = f° k , g k = go9 k . 



(fk , E k+ i(A)A k+ ig k ) dfi x dfi(x) = / / (/* , AA k+1 g k ) dfi x dfi(x) 
: (f k , A k+ ig) = (f k , P k+ iAP k+ igk) = (fk , Ag k ) = (fk , A k gk) 



(f k , AA k g k ) dfi x dn(x) = / (f k , E k (A)A k g k ) dfi x dfj,(x). 
xJsi, JxJn x 

So for /x-a.e. x e X, E k {E k +i{A)A k+ i)(x,-) = (E k (A)A k ){x,-), fi x -a.e., for all k > 1. Therefore the 

sequence {(.Efc(A).Afc)(x, •)}/. is a martingale. By Doob's martingale convergence theorem, Ek(A)Ak(x, ■) 

converges pointwise fj, x -&.e. to AA(x, •). 

Now we compute Ek(A)(x, •). We have 

(5-2) / (fk , E k {A)g k ) dfi x = [ (f k , Ag k ) dfi x = 1 V (f(y) , m< fc )*fc(y)m<*> (y)g(y)) . 
Ja x Jn x #r k {x) £-> \ I 

Let C Uli ... iUJk be the cylinder of points in fl x that start with LOi,....uj k . Let f k and g k be supported on 
Ctt)i,...,ai fc - Then we obtain from ()5.2() that 

£ , fc(A)(x,Wi, . . . ,Wfc)Ar(C'wi,...,w*) = # r -k( x ) (™ (fc) *^™ fc ) (Wfe)- 

But, with Lemma [5~4l 



£x(C Ul ,...,«J = Tr(P x (CU,..., w J) - Tr f - _\ (mW'kW)^) 

(m^*hm^) o 6> fe 



Thus 

(5.3) Bjk(A) = 

Then 



E k (A)A 



k 



Tr((mW /imW) o 



Tr((m( fc ) /im( fc ))o6» fe ) 

This proves the theorem. □ 
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Remark 5.7. In the scalar case, the terms involving my 1 ' will disappear and we reobtain the results from 
|DJ06bj . 

Remark 5.8. An ergodic limit for low-pass filters. Our limit theorem can be used to obtain an 
interesting limit result. The matrices enable us to use a trick to compare two different measures. 

Let X = T 1 , r(z) = z 2 . Denote by ipi,ip 2 the inverse branches of r, ipi{e ie ) = e ie / 2 , ipi(e ie ) = e l ( e+27r )/ 2 , 
for 9 £ [0, 2tt). Then the solenoid X^ is measurably isomorphic to X x Q, where f2 := {1,2} N , because a 
point (zq, zi, . ■ .) € Xoa consists of zq G T 1 and a choice of the inverse branches u>x, u>2, ■ ■ ■ ■ 

Let mi(z) = (1 + z)/v2 (or any low-pass filter that gives orthogonal scaling functions in L 2 (M)), and 
mn{z) = 1. Let h% = hi = 1. Then the covariant representation for (mi, hi) is P 2 (R) so the measure P x 
on the solenoid is supported on sequences u € such that there exists uq such that u) n — for all n > no 
or bj n — 1 for all n > no ( scc |DJ06a , or Section [5]). The covariant representation for (m2,/i2) is on the 
solenoid X^ with the Haar measure and the measures P 2 are the Bernoulli measures on Q, where {0} 
and {1} have equal probabilities 1/2 (see also [DJ06a and 
mi 



Let m 



m 2 



Then the measure P x is clearly 



DJ06cQ . 

" Pi 



P 2 



i.e., for / e C(X Q 



fdP x = 



f Xoo fdPid»(x) 




X. 





fdP 2 d^(x) 



The trace of this itigrs lire is jl x — -\- P^ . As we explained before (see also Section the measure P x 
is atomic, and let us denote the support of this measure by R x . On the other hand, P 2 is the Bernoulli 
measure, so it has no atoms, therefore R x has P 2 -measure zero. 



Then it is easy to see that -j^- = xr. 



A 



and -p^- — 

dP x 
dfi x 



Xn\R x 

XR* 





Then 




Note that this shows that A (a;) is singular everywhere. 

We will remark that this implies an interesting phenomenon, which occurs for any low-pass filter that 
gives orthogonal scaling functions in L 2 (R). 

First, as shown in |D J06a] . since P x is supported in on R x , it follows that for any to outside R x one has 



lim 



777, 



However, from the convergence Theorem I5.6[ we have that for /x-a.e 

mW*mW(^,...(/) Ul i) 



lim 

k — ^oo Tr(7Tl' ' 

This implies that for fl x a.e. u> € Q \ R x 

= lim 



m ( - kS >(ip ulk ...ip LUl x)) 



\m[ k \ip UJk ...^ UJl x)\ 2 



= Pi(M) = o. 

x € T 1 and for fl x -a.e. to: 
A(x,uj). 



in 



(fe)/ 



so limfe^oo 77i (fc) (i/} Uh . 



o 1 (ip u)k ..4 UJl x)\ 2 + 1 

.ip UJl x) = for fl x -a.e. w, so we obtain the much stronger limit for P 2 -a.e. ui (and 



recall that P x is the Bernoulli measure): 



lim m i - k \ip UJk ...^ LJl x) 



0. 



6. Low-pass filters 

In the scalar case, if the filter tuq satisfies a low- pass condition mo(l) = -\/2, the classical wavelet theory 
shows that the scaling equation has a solution in L 2 (R). Of course the solution might be a non-orthogonal 
scaling function, and then there are super- wavelet constructions (see [BDP05 ) that will give orthogonal 
solutions, functions. 

The point we want to make is that, when a low-pass condition is satisfied, the resulting covariant repre- 
sentation can be realized on L 2 (R), as in classical wavelet theory, or in a direct sum of copies of L 2 (M), as 
in the super- wavelet theory developed in BD P05] . 
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In the matricial case, the low-pass condition is replaced by the -E(7)-condition introduced in [JS99J . 

Our covariant representations are on the solenoid A^, but we show that when a low-pass condition is 
satisfied the measures P x are supported on an embedding of R in the solenoid, and P x is directly related to 
the scaling functions. 

We recall now the setup from [DR06] . We assume that X is a complete metric space with an isometric 
covering space group action of a group G such that X/G —: X is compact. Moreover, we assume that 
f : X — > X is a strictly expansive homeomorphism and there exists an endomorphism A G End(G) such that 
fg = (Ag)f and AG is a normal subgroup of index q. Since f is expansive it has a fixed point xq. 

Let p : X — > X denote the quotient covering map and define r : X — ► X as r(p(x)) — p(fx). Let 
x o '■= p{xq). Let fi be a strongly invariant measure on X and fi the measure on X obtained by lifting the 
measure \i by the covering map p (see also (16. 1ft ) . 

We assume that the bundle £ over X is a Lipschitz continuous bundle, and that p*£ is trivial. 

Example 6.1. The main example is the one used in wavelet theory: X = R™. The group G = Z" acts 
on 1" by translations: g, x i— > x + g, (x G R", g G Z n ). Define the map f(x) — Ax, where A is an n x n 
expansive integer matrix. The quotient R™/Z™ can be identified with the torus T™, p : R™ — > R n /Z n is the 
quotient map, and let r(x) = Ax mod Z™ for x £ T". The fixed point of f is xo = 0, and Xo := p(xo) = 0. 
The bundle £ over T" is T" x C d , and p*^ = R™ x C d . 

We let m € Af ( ;(Lip 1 (X)). Moreover, we assume that 

• Rl = 1; 

• The matrix m(xo)/ v / g satisfies the -E(7) condition (according to JS99J, i.e., 1 is the only eigen- 
value of m{xo) / yfq with absolute value greater than or equal to 1, and its algebraic and geometric 
multiplicity are both equal to I > 1. Let i?i denote the eigenspace corresponding to the eigenvalue 
1. 

Recall that S is the set of continuous sections in the bundle £. 

Let S := {/ e C h (X)| J] geG |/| 2 ° 5 e C(X)}. S is a C(X)-Hilbert module with the inner-product: 

This module inner- product is related to the ordinary inner-product from L 2 (X,j2), by 
(6.1) f (r,dfi= f (CvYdfx, (CtjGS). 

Let [/ G B(L 2 (X,fl)) denote the unitary operator defined by 

Uf = q 1/2 f o f, and , Ui := U@ ■ ■ ■<&& . 

i times 

Define the representation tti : C(X) — > B(® l j =1 L 2 (X , fi)) by 

7f,(o)(/i, . . . , /,)(x) = {a(px)h{x), a(jpx)fi(x)), (a G C(X)J U ...,/; G L 2 (l,/i)). 

We will need to define some "scaling functions". These will be fixed points of a refinement operator 
obtained as the limits of the iterates of this refinement operator. 
To initialize the iteration, we fix an / G S such that: 

• </,/)' = ! 

• (Z g£G \f(9x)-f(gy)\ 2 ) 1/2 <Dd(x,y) 

• f(gxo) = for every g ^ 1. 

Let si, . . . , si G S be Lipschitz sections in £ such that Si(xo), . . . , s/(xo) form an orthonormal basis for E±. 
As in |DR06| . we define the starting points for the cascade algorithm Wj G Hom C ( X ) (S, S) by W}s = 
(sj, s)f and the refinement operator 

M : Komc( X )(S, E) -> Hom c(x) (5, H), (Mf)s = r¥ms o r, for G Hom c(x) (,S', H), s G 5. 

Proposition 6.2. [DR06] 
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(i) The following limit exists and is uniform on compact sets: 

V{x) := lim q- k/2 m ( ' k) (p(r- k x)), (x G X). 

k — >oo 

Also, V(xq) is the projection onto the eigenspace E\. 

(ii) For each j S {1, . . . ,1} and for s£ S, {M k Wjs}k>i converges uniformly on compact sets to WjS, 
and this defines Wj £ Homc<x){S^'S), Wjs(x) = (sj(xo) , V(x)s(px)} , s G S , x £ X . Moreover 
MWj = Wj. 

(iii) For each j S {1, . . . , 1} the map hj :— W*Wj defines a minimal projection in the algebra of contin- 
uous harmonic functions Sj c (hj is a projection also in the algebra of bounded measurable harmonic 
functions 9) but it is not necessarily minimal). 

hj(px)s(px) = ^2 (P*(9 x ) s j( x o) , s(px))'P*(gx)s j (xo), (s £ S,x G X). 
g&G 

Moreover these projections are mutually orthogonal in this algebra. 

Theorem 6.3. Let h :— hj an d ^ % denote Hilbert space of the covariant representation obtained 

from m and h. There exists a unitary J : Tt — ► (Bj=i L 2 {X,fi) such that 

(i) UiJ = JU 

(ii) Js o 6>o = WjS for every s £ S. 

(iii) iti(a)J — Jir(a) for every a £ C(X). 

Proof. Define Jk ■ Tit —> © 3 L 2 (X, [i) by JkU~ h s o 9 = © ■ U~ k WjS. Jk is an isometry because 

(J k U- k so6 Q , J k U- k so6 ) L 2 =^(?7- fc W 3S ,C/- fc W, S ) L 2 = ^{W j s,W j s) L ^=Y J f (W jS , V^s/d/i 

= I ( s >Yl = / <* , hs) d/i = (so O , s o 9 Q ) = (U- k s o 9 , lJ- k s o Q ) 

JX j Jx 

Moreover, since Wj — MWj — U~ x Wjms o r, 

JkU-^-^s o9 a = J k U- k Us o 9 = ® j U- k W j ms or = ej£/" -( * _x) VV,-s = J fc _iZ7 _(k-1) s o 9 

so Jk\H k -i = Jk—i- This gives us an isometric map J : Ji — > j L 2 (X , ft) , and the intertwining properties 
of J are checked by a direct computation. 

Let H := UkJH-k C © 3 L 2 (X,fi) and let Q denote the orthogonal projection onto this space. We have 
UiQ = QUi, Tii{a)Q = Qfrj(a) for every a € C(X). 

We have that Uf k TT l (a)U l k is a multiplication by the matrix that has aopo f k on the diagonal. Since f 
is expansive, an application of the Stone- Weierstrass theorem shows that {aopo r~ k | a £ C(-X"), fc > 0} is 
dense in C C (X). Since Q commutes with all operators of the form Uf k ni(a)U k , this implies that Q commutes 
with L°°(X,p,) (8> Since L°°(X,fi) is a maximal abelian subalgebra, it follows that Q corresponds to a 
pointwise multiplication by a map in L°°(X,ji) ® M;(C). Since Q is a projection, this map is projection 
valued. Moreover, Q(x) — Q(rx) fi-a.e., since QUi — UiQ. 

Consider now ifi :— ®™ =1 VVjSj for i £ {1,...,/}. We have with Proposition 16.21 <fi(xo) = e», the 
canonical vectors in C'. Then, using the continuity of Wj, we have that for i in a neighborhood of xq, 
{(Pi(x) \ i £ {1, . . . , I}} forms a basis for C l . Since ipi £ H, Qipi — cpi so Q(x)<pi(x) = tpi(x) for all x £ X. 
But then, Q(x) must be the identity in a neighborhood of xq, and since Q(x) = Q{rx), and f is expansive, 
we obtain that Q(x) is the identity for all x € X. Thus Q = 1 and is the entire space ® ■ i 2 (X, /i). 

□ 

Proposition 6.4. T/ie map i : X — > Xoo, i(x) = (p(f _ a;))^.]^ is a continuous bijection onto the set of 
sequences (Zk)k G ^oo wit/i limfc^oo = xq. For all p(x) £ X , i(X) D Q. px = i(Gx). 

Let P x be the measures associated to h — YZj hj as in Proposition \5.3\ Then P px is atomic and supported 
on i{Gx), and 

P px ({l(gx)}) = V{gx)*V{gx), (x £ X). 
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Proof. Since r is expansive, the sequence f~ k x converges to the fixed point xq, so for k large r~ k x is in some 
neighborhood where p is injective. This implies that i is injective. 

The continuity of i is clear, and p(f~ x) converges to p(xq) = Xq. To see that i is onto the given set, 
take some sequence (zk)k in Xoo such that Zk converges to xq- Take a neighborhood V of xq such that the 
restriction of p to V is a homeomorphism onto the neighborhood p(V) of xo- Take a smaller neighborhood 
U C V of Xq such that f(U) C V. For fc large z k is in p(J7). So z% — p(xk) for some Xk £ U. Since 
f{zk+i) — Zk it follows that f(xk+i) — gxk for some g £ G. But as both f(x/c+i) and X& are in V, it follows 
that g must be 1. So Xk.+i = r~ 1 Xk. for k large, bigger than some ko. Then if define x :— r~ k °Xk , we have 
i(x) = {z k )k- 

If i(y) is in fl px then py — px so y — gx for some g £ G. 
Let us check the P^-measure of the atoms. We have 

Ppx{{i{gx)}) = lim P pa; ({(z„) rl | Zj =p(r _J 'x),0 < j < fc}) = lim q- k m {k) * {pf- k x)h{pf- k x)m {k) {pf- k x) 

k— *oo A;— >oo 

= P(x)*/i(x )P(x) = P(x)*P(x), 
and we used the fact from [DR06] that the range of ^(x) is contained in E\ and /i(xo) = Y^j hj( x o) ^ s the 
projection onto E\. 

Take the cylinder C px<zl> ,,, tZn of sequences in r2 p:c that start with px, zi, . . . , z„. If we add all the atoms 
in this cylinder, we obtain, 

(6.2) Yl P P Mgx))= Y, T(x + g yV(x + g). 

g,l{gx)eC px ,g! 1: ...z n S>i(sa:)eC'p a ,, iS1 ,... an 

On the other hand , for any section s e S, and with the notation px n := z n , for some x n G X, using the 
formula for h in Proposition 16. 21 

(s(z ), P P x(Cp X . Zu ..., z Js(z )} = (s(z ), q~ n m {n) *{z n )h{z n )m {n) {z n )s(z fj 

= J2T,^ n (m [n) (z n )s(z ) , (v*(gx n ) Sj (x ) , m (n \z n ) s (z )) V*(gx n ) Sj (x )) 
j gee 

= Y,T,\( V *^x n )s 3 (x ) , q-^ 2 m^(z n ) S (zo))\\ 
g j 

and since {sj(xo)}j is an orthonormal basis for E\ and the range of V(x) is contained in E\, 

= Y,\\ V (9x n )q- n/2 rn (n) (z n ) S (z )\\ 2 . 
g 

But, from the definition of V, we have T'{gx n )q~ n ^ 2 m! yn \z n ) = V(g'x) for some g' € G with i(g'x) £ 
Cp X ,zx,—,z„i and we obtain further 

\\V{g'x) s {z )f. 

Comparing with (|6.2[) . this shows that the sum of the atoms is equal to the measure of the cylinder and 
thus the measures P px are supported on these atoms. □ 

Remark 6.5. Scaling functions. To obtain old-fashioned scaling functions as described in the introduc- 
tion, let us consider the case when we are dealing with the Example 1 6. 11 used in the regular wavelet theory 
(but the arguments below work also in the more general case we described in this section) , and let us consider 
the case when £ is the trivial vector bundle X x C d = T" x <C d . 

Then we can take the canonical sections in £, Cj(x) =e„ieT", ie{l,...,d}, where ej are the canonical 
vectors in C d . Then, define for each j £ {1, . . . , I}, 

<fi:=WjCi, (i £ {l,...,d}). 

Then ipl is in 5 so it is a function in L 2 (K d ) (according to (|6.1|1 ). Also, since MWj = Wj, we have 

d d 

U(pl = UMWjCi = WjiriCi o f — WjQj m ki c k ) = 2J H m hi) l fi- 

k=l k=l 
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(Here n(s)f(x) = s(px)f(x), s G C*(T"), / G L 2 (R"), x G 1"). 
Also, for / G C(T n ), 



(vl , <SWr) = I ¥j4 dx= f f U , 4)' dn= f f (W jCi , W jCi ,) dn= f f (a , W*W jCi ,) dn 



/ (cj , hjCi>) dfi = f(hj)u>dfi. 

Thus, ip{, . . . , ip J n form a multi-scaling function in L 2 (R"), with filter m and correlation matrix /ij. 
We can put together all these multi-scaling functions and define 

<Pi := (<pl . . . , <4) G ®< =1 L 2 (K"), (i G {1, . . . , d}). 

Then we still have the same scaling equation, but now in ©^ =1 L 2 (R™): 

d 

Uifi = /,n{mki)<Pk, (i G {1, . . . , d}), 
fe=i 

and the correlation matrix for tp±, . . . , tpd is the harmonic function /i, i.e., 

(<Pi,n(f)<Pi>) = I fhi'dfi, (/gC(TT)). 

In Theorem 16. 31 we see that if s — (s 1 , . . . , s d ), then 

d d d 

j S oe = (B^Wjs = eJ =1 Wj(E s<c i) = ©!=i(E^( sl )^) = E^( fil )^- 

i — 1 i—1 2—1 

This implies, with Theorem 16.31 that the linear span of 

{Ufa(f)<Pk 1 3 e Z, / G C(T"), fc G {1, . . . , d}} 

is dense in ® l j=1 L 2 (R n ). 

Thus, ipi, . . . ,ipd form a "super" -multi-scaling function for the bigger- ( "super-" )space ®j =1 L 2 (K"). 

Acknowledgements. We would like to thank professor Palle Jorgensen for his suggestions and for offering us 
a wider view on the subject, which we included in the introduction. 
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